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ABSTRACT 

The CLASH X-ray selected sample of 20 galaxy clusters contains ten brightest cluster galaxies (BCGs) 
that exhibit significant (>5 a) extinction-corrected star formation rates (SFRs). Star formation 
activity is inferred from photometric estimates of UV and H<a+[N II] emission in knots and filaments 
detected in CLASH HST ACS and WFC3 observations. UV-derived SFRs in these BCGs span two 
orders of magnitude, including two with a SFR > 100 M 0 yr _1 . These measurements are supplemented 
with [O II], [O III], and H/3 fluxes measured from spectra obtained with the SOAR telescope. We 
confirm that photoionization from ongoing star formation powers the line emission nebulae in these 
BCGs, although in many BCGs there is also evidence of a LINER-like contribution to the line emission. 
Coupling these data with Chandra X-ray measurements, we infer that the star formation occurs 
exclusively in low-entropy cluster cores and exhibits a correlation with gas properties related to cooling. 
We also perform an in-depth study of the starburst history of the BCG in the cluster RXJ1532.9+3021, 
and create 2D maps of stellar properties on scales down to ^350 pc. These maps reveal evidence for 
an ongoing burst occurring in elongated filaments, generally on ~ 0.5-1.0 Gyr timescales, although 
some filaments are consistent with much younger (< 100 Myr) burst timescales and may be correlated 
with recent activity from the AGN. The relationship between BCG SFRs and the surrounding ICM 
gas properties provide new support for the process of feedback-regulated cooling in galaxy clusters 
and is consistent with recent theoretical predictions. 

Keywords: galaxies: clusters: general - galaxies: clusters: intracluster medium - galaxies: starburst 


1. INTRODUCTION 


Brightest cluster galaxies (BCGs) in cool core galaxy 
clusters exhibit nebular emission features that are 
thought to be related to the intraclu ster medium ( ICM) 
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some of this activity is triggered by active galactic nu¬ 
clei (AGN), most of the emission appears to be powered 
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been observed m the Phoenix cluster (z = 0.596), which 
hosts a massive starbursting BCG producing ne w stars 
at a r ate of nearly 1000 M e yr -1 ( McDonald et al.|2012 


2013). 

The most plausible candidate for the source of the gas 
being converted into stars in these galaxies is radiatively 
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images of low-redshij 

ft (z^O.l) BCGs suggest that while 


cooled ICM plasma (Fabian|1994|. In a cool core cluster, 
the cooling time below a critical radius is less than the 
Hubble time at the redshift of the cluster, and plasma ini¬ 
tially at this radius ought to have cooled. The ICM that 
manages to cool descends into the gravitational well in 
order to maintain pressure equilibrium, ultimately con¬ 
densing in to star forming gas inside the virial radius of 
the BCG (Fabian 19941). However, the cooling inferred 
no^v 


the s tar formation rates i n these syste ms (e.g. 
et al.|1989* 


from a simple cooling How model is far more rapid than 
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Activity m BCGs in the form of recent star formation 
and AGN outbursts are important components of feed¬ 
back mechanisms proposed to reconcile the tension be¬ 
tween the pre dicted and o bserved ICM coo ling in c_ool 


2015]). Specifically, the hot, X-ray emitting ICM of a re- 
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Taxed galaxy cluster is predicted to radiatively cool more 
rapidly than is typically observed, and this manifests in 
star formation rates (SFRs) in BCGs which are roughly 
an order of magnitude lower than what would be pre¬ 
dicted if all the available ga s did indeed condense into 


cold gas (O’Dea et al. 2008). However, in the presence 


of a feedback mechanism, such as energy injection from 
an AGN, the ICM is partially reheated and prevented 
from cooling catastrophically. Instead, residual cooling 
or cooling during an off-mode in the feedback duty cycle 
can account for the extended star forming structures ob¬ 
served. Studying the residual cooling using observations 
of BCGs spanning a wide range of activity will allow 
us to learn about different phases of cooling and feed¬ 
back, and help us to determine whether a single feedback 
mechanism accounts for the variety of BCG features we 
observe. 
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The high-resolution, multi-band HST observations 



galaxy clusters, and found evidence for significant, ex¬ 
tended emission attributable to recent star formation in 
10 of them. Two BCGs in this sample, RXJ1532.9+3021 
and MACS1931.8-2653, stand out due to their strikingly 
large and luminous UV filaments. 

In the present study, we examine BCGs in the subset of 
CLASH clusters that were X-ray selected. This subsam¬ 
ple includes 20 of the 25 CLASH clusters and includes all 
of the star-forming BCGs as identified by UV features. 
Using the CLASH HST photometry, along with spec¬ 
tra from the Southern Astrophysical Research (SOAR) 
telescope and archival Chandra data, we investigate the 
nature of star formation in these BCGs, and provide new 
constraints on the source of the star formation activity 
in the structures we observe. 

This paper consists of two parts. First we derive 
SFRs for all BCGs using HST photometry and charac¬ 
terize the source of nebular line-emission in UV-bright 
BCGs using a combination of the HST photometry 
and SOAR-Goodman spectra. Second, we analyse the 
connection between the star formation and the prop¬ 
erties of the ICM. This second part includes a de¬ 
tailed star formation history (SFH) analysis of the BCG 
in RXJ1532.9+3021 derived from SED fitting of the 
CLASH photometry to create maps of stellar population 
parameters. RXJ1532.9+3021 was chosen for more de¬ 
tailed study because of the spectacular nature of its UV 
and Ha structure (see Figure [l]). 

We report on the incidence and distribution of 
reddening-corrected SFRs in our sample, and demon¬ 
strate that structures qualitatively comparable with 
those in the massive outburst in the Phoenix cluster are 
the sites of BCG activity in these ‘intermediate’ star- 
bursts as well. Structures observed in CLASH clusters 


also be ar similarities to BCGs analysed i n [Tremblay et al. 
(2015). These SFRs are compared with Chandra derived 
KJ1V1 core entropies and predicted cooling rates in the 
low-entropy core ICM in order to test the hypothesis that 
star formation is being fed by ICM cooling. 

For RXJ1532.9+3021 we are able to compare the prop¬ 
erties of BCG filaments with the properties of X-ray cav¬ 
ities in the ICM, which we use to assess recent AGN 
jet-mode activity. We compare our results to pr evio us 


analysis of the morph ology of the cluster ICM (Hlavacek- 
Larrondo et al.|2013 ). The maps produced for this BCG 


allow us to examine the star formation history in indi¬ 
vidual knots and filaments, down to ~ 350 pc scales. 
These data allow us to investigate the source of the gas 
that condenses into the star forming regions and, for the 
case of RXJ1532.9+3021, to examine this condensation 
in significant detail. 

Our paper is organized as follows. We describe the ob¬ 
servational data in Section 2. In Section 3, we present 
the analysis of these data. Results are presented in Sec¬ 
tion 4, and the astrophysical implications are discussed 
in Section 5. We summarize our main conclusions in 
Section 6. We adopt the following cosmological param¬ 
eters throughout this work: H 0 = 70.0 km s -1 Mpc -1 , 


= 0.30, and = 0.70. We assume a Salpeter (1955) 
IMF throughout. 


2. OBSERVATIONS 
2.1. CLASH HST Observations 


The CLASH program is detailed in Postm an et al.| 
(2012J. The 20 X-ray selected clusters were each ob¬ 
served for 15-20 orbits (for a total of at least 20 orbits 
including archival data) divided among 16 bands of pho¬ 
tometry spanning an observer-frame wavelength range of 
^2000-17000 A. We use multi-band mo saics drizzled to 
a common 07065 pixel scale (Koekemoer et al.| 2011). All 
flux measure ments are corrected fo r foreground redden- 
ing using the (Schlegel et al. 1998) dust maps. A single 
reddening correction due to dust m the Milky Way was 
calculated for each BCG in each filter. Additionally, we 
perform a background subtraction using a combination 
of iterative 3-sigma clipping on large scales along with 
a local median flux measurement in an annulus around 
each BCG in the UV. 

Early CLASH WFC3/UVIS observations are affected 
by non-uniform flat-fielding on scales of hundreds to 
thousands of pixels. We accounted for this in observa¬ 
tions of BCGs with faint or possibly no significant detec¬ 
tion of UV by extracting photometry from multiple iden¬ 
tical apertures in WFC3/UVIS images placed on empty 
patches of sky, and adding the scatter in these apertures 
to our error budget. For all CLASH observations, we 
apply a 3% floor to the total photometric uncertainty to 
account for all sources of systematic error and absolute 
flux calibration uncertainty. 

2.2. Chandra X-ray Observations 

All the clusters in the CLASH X-ray sample have 
archival Chandra data. For this work, we use the ICM 
parameters published in the Archive o f Chandra Clus¬ 
ter E ntropy Profile Tables (ACCEPT) (|Cavagnolo et al. 
20091). Gas density and temperature profiles were mea¬ 


sured in ACCEPT, and the core entropies and cooling 
times used in the present study are calculated using these 
profiles. Gas density profiles were measured in concentric 
annuli 5” wide. Temperature profiles, which were used 
in combination with the gas density profiles to derive 
both cooling time and entropy profiles, were measured 
in concentric annuli containing at least 2500 counts per 
bin. 


2.3. Spectra 

Optical spectra of the BCGs were obtained for 15 of 
the 20 CLASH X-ray selected clusters. Objects were 
observed with the Goodman High Throughput Spectro¬ 
graph on the SOAR 4.1 meter telescope using either the 
KOSI 600 grating or the SYZY 400 grating. The KOSI 
grating’s dispersion is approximately 0.65 A pix -1 , while 
the SYZY grating’s is approximately 1.0 A pix -1 . Their 
spectral ranges are roughly 2670 A and 4000 A, respec¬ 
tively (see Table 1). Central wavelengths were selected 
to best include the [O II] doublet ([O II] A,A3726,29), 
the [O III] doublet (O III] A, A4959,5007), and H/3 at 
their redshifted positions. Position angles were chosen to 
sample observed filamentary structures or other objects 
of interest. For all observations, the 1.68" long slit was 
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Figure 1 . Color composite images of the most active star-forming BCGs of RXJ1532.9+3021 (z=0.363) on the 
left and MACS1931.8-2653 (z=0.352) on the right. The RGB color composites are made using WFC3-IR filters 
F105W+F110W+F125W+F140W+F160W in red, the ACS filters F606W + F625W + F775W + F814W + F850LP in green, 
and the ACS filters F435W + F475W in blue. 


used. 


Obs ervations were reduced following | Werner et al. 
(2014). Bias correction and trimming were performed 
with the IRAfljtask CCDPROC. Quartz lamp frames, taken 
before and after observation images, were used to flat- 
field the images. Wavelength calibration was performed 
with FeAr arc lamp exposures, with distortion along the 
spatial direction corrected for by tracing the position 
of standard stars. Sensitivity functions were produced 
from same-night observations of standard stars with the 
10.0" long slit that, along with an extinction correction 
based on IRAF’s extinction file ctioextinct.dat and 
a correction for air mass, were used to flux calibrate 
observations. After background subtraction to minimize 
the contribution of night sky lines, observations were 
then median combined using IMCOMBINE. 


3. ANALYSIS 

3.1. Mean UV Luminosities 
We calculate the UV lumino sity of C LASH BCGs in or 


der to estimate SFRs using the Kennicutt (1998) relation. 
To do this, we extract the flux from those CLASH filters 
whose pivot wavelengths fall in the range 1500 —2800A in 
the cluster rest frame , correspo nding to the wavelength 
range used by Kennicutt (1998) to calibrate the L uv~ 
SFR relation. Si nce the continuu m L uv due to young 
stars is flat for a Salpeter (1955) IMF and continuous 
star formation, we use fluxes extracted from the relevant 


1 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation. 


UV filters to calculate the average luminosity (L v ^jjy). 
This quantity is one of two SFR proxies we are able to 
calculate from the HST photometry alone, the other be¬ 
ing L hol+[nii]i which we discuss in Section 3.2. 

For the majority of the BCGs, we extract fluxes from 
a subset of the CLASH WFC3/UVIS filters, F225W, 
F275W, F336W, and F390W. For the highest redshift 
clusters (z > 0.542), we use ACS filters F435W and 
sometimes F475W as well. Filter selections are shown 
in Table 2. 

Since we are interested in the UV flux from young 
stars, we must remove the contribution from the UV- 
upturn in the quiescent stellar population in the BCGs. 
The UV-upturn has been well-studied and is due to post 
AGB stars an d e xtreme blue horizontal branch stars (e+ 


Brown| |20041 |Yi||2008| |Ferguson fc Davidsen||1993l |YI| 
et al. |1998| ). T his component of the UV flux mimics 


low-level star formation and would, if left uncorrected, 
bias our SFR estimates in UV-faint BCGs. 

A template SED of the UV-upturn population is de¬ 
rived by averaging the normalized ACS-WFC3IR rest- 
frame photometry of 5 satellite galaxies in each of 20 
CLASH clusters (for a total of 100 satellite galaxies). 
We do not include WFC3/UVIS in the satellite galaxy 
photometry to avoid complications in estimating the sen¬ 
sitivity to weak UV sources, and instead opt to use the 
GALEX-2MASS J colors in Hicks et al. (2010) to extend 
the range of our data below 2000 A rest . These data are 
fit to a fifth order spline, and the uncertainty in the fit is 
estimated using a Monte Carlo distribution of spline fits. 

We estimate the UV contribution from old stars by 
subtracting out a model of the underlying early-type 
galaxy. A 2D non-parametric model of each BCG was fit 
to the F160W image, which is dominated by the BCG’s 
old stellar population. The model is scaled according to 
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Table 1 

SOAR Observations of Brightest Cluster Galaxies 


Cluster 


Obs Date Exposure Times Grating Range PA a Airmass Standard Star 

(YYYY-MM-DD) (s) (1/mm) A ( °) 


Abell 209 

2012-09-24 

1x1200, 1x1800 

400 4514-7555 

135 

1.2 

LTT7379 

Abell 383 

2012-10-09 

2x1200, 1x600 

400 4126-7568 

0 

1.3 

LTT1020 

MACS0329.7-0211 

2012-11-09 

3x1200 

400 4612-7555 

125 

1.1 

LTT7379 

MACS0429.6—0253 

2012-11-19 

2x1200, 1x900 

400 4566-7556 

167 

1.3 

LTT1020 

MACS1115.9+0219 

2013-05-11 

4x900 

400 4401-8462 

130 

1.2 

LTT4364 

MACS1206—0847 

2013-05-11 

4x900 

400 4398-8457 

100 

1.1 

LTT4364 

MACS1311.0—0310 

2013-05-11 

1x1500, 1x900 

400 4412-8470 

40 

1.5 

LTT4364 

MACS1423.8+2404 

2015-02-26 

3x1200 

600 5049-7724 

0 

1.8 

LTT4364 

MACS1720.3+3536 

2015-06-13 

3x1200 

600 4570-7235 

160 

2.2 

LTT6248 

MACS1931.8-2653 

2012-04-17 

1x1800 

600 4271-6938 

252 

1.0 

LTT7379 

MS2137—2353 

2015-06-13 

4x1200 

600 4570-7235 

145 

1.1 

LTT6248 

RXJ1347.5—1145 

2012-07-15 

1x1200, 1x900 

600 4556-7223 

125 

1.6 

LTT9491 

RXJ1532.9+3021 

2012-04-17 

2x1200 

600 4811-7471 

187 

2.1 

LTT3864 

RXJ2129.7+0005 

2012-07-15 

4x1200 

600 4560-7228 

201 

1.2 

LTT3864 

RXJ2248.7-4431 

2013-09-08 

4x1200 

600 4299-6969 

352 

1.1 

LTT1020 


a Position Angle measured east of north. 



Figure 2. Rest frame SED for quiescent starlight in BCGs. 
The blue points are ACS/WFC3IR photometry taken from 
100 satellite galaxies, and represent the averages of 5 galaxies 
per CLASH cluster in each of 20 clusters. The red points are 
derived using the GALEX-2MASS J color calibration from 
Hicks et al. 2010 for a sample of quiescent BCGs at an average 
redshift of 0.115. The FUV/NUV rest wavelengths depicted 
in this figure are for Abell 209 (z m 0.209), after accounting 
for the average redshift of the Hicks et al. color calibration. 
The grey line represents the spline fit to the data. Fluxes are 
reported in normalized units of F u . 

the template SED shown in Figure[2]and subtracted from 
the UV images. The residual UV flux measured after the 
subtraction is our UV-upturn-corrected estimate due to 
ongoing star formation activity. 

3.2. Ha+[N II] Maps 

The CLASH data allow us to estimate Ha emission in 
the BCGs using a broadband subtraction technique in 


cases where the emission is strong relative to the stellar 
continuum. For each BCG, a ‘line’ and ‘continuum’ filter 
is chosen based on the cluster redshift (see Table 2). One 
or more satellite early-type galaxies are chosen based on 
having an IR—optical color similar to the BCG, and the 
mean ratio of the satellite galaxy fluxes between the line 
and continuum filters is used to scale the continuum fil¬ 
ter to the line filter. We subtract the scaled continuum 
filter image from the line filter image, leaving a residual 
flux that is primarily due to Hct+[N II] emission. The 
Ha nebulae of clusters with significant UV emission are 
shown in Figure [3J 

This method for estimating iLct+[N II] is reliable only 
for line emission with a large equivalent width (EW). 
iLct+[N II] EWs may be approximated by 


7 \ -o /line /cont /-,\ 

dA^Bii ne - - -, (1) 

Jcont 


where f\[ ne and / CO nt are the fluxes through the line and 
continuum filters, and Bii ne is the photometric band¬ 
width of the line filter. If we assume a 3% uncertainty for 
the surface photometry in both the ‘line’ and ‘continuum’ 
filters, we can reliably recover features with iLct+[N II] 
EWs that are 0.044 xBi ine . Because the BCGs are very 
bright in both the ‘line’ and ‘continuum’ filters, so the 
dominant source of uncertainty on our photometry will 
be the absolute flux calibration. For BCGs where we 
adopt ACS filters (either F775W or F850LP) as the line 
filters, we can recover iLct+[N II] features with EWs > 
20 A in the observer frame. For BCGs where we adopt 
WFC3IR filters (either F105W or F125W), we can re¬ 
cover features with EWs > 40 A in the observer frame. 


EW = 


/ 


JX /cont 
/cont 
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3.3. Other Emission Lines 

Longslit spectra provide coverage of the [O II], 
[O III] A, A4959, 5007, and H [3 emission lines. We mea¬ 
sure line luminosities by fitting a Gaussian line profile 
and a continuum to reduced, 1-D longslit spectra using 
the IRAF task splot. Continuum levels in splot are 
identified by averaging regions of continuum emission 
adjacent to emission lines, and continuum-subtracted 
Gaussian line profiles are fit using the default iterative 
Levenberg-Marquardt algorithm. An [O II] luminosity 
could not be calculated for Abell 383, owing to contam¬ 
ination of the spectrum at ^ 4424 A. [O III] A 5007 is 
unavailable for RXJ1347.5-1145 since the available spec¬ 
trum does include this line. 

[O II] luminosities are an independent check on the 
SFR derived from our UV and Ha luminosities. While 
[O II] luminosities are usually considered to be a less reli¬ 
able estimator of SFRs due to their dependence on ioniza¬ 
tion and metallicity, they can provide extra constraining 
power when used in con junction with other SFR est i¬ 
mators ([Chariot fc Longhetti 2001; Rosa-Gonzalez et al 


2002 Kewley et al 


s Long. 

.2004; 


Moiistakas et al 


We can 

use luminosities in "a similar fashion. Assuming case 
B recombination, we can de rive an S FR fr om H/?, using 
the relation Ha/H/3 = 2.85 (fVeilleux 2002). 

[O III] A 5007 is useful for constructing diagnostic dia¬ 
grams to separate regions heated by normal stars versus 
other sources of ionization like AGN and shocks. The 
classic diagram for distinguishing star forming regions 
from AGNs is the BPT diagram ^Baldwin et al. 1981 


Kewley et al. 


2001a Kauffmann et 


TTpOOdbi). Wecon- 
struct a modified version of the BPT diagram based on 
a combination of spectral and broadband data. We use 
an additional diagnostic, when possible, that is nearly 
insensitive to extinction, and only uses our spectroscopic 
data. The so-called ‘blue-line’ diagram barely depends 
on the accuracy of the reddening correction because it 
is derived solely from equivalent width values. In this 
work, we use the ‘blue-line ’ diagram derived from the 
[O II], fO II I], and H/3 lines (|Lamareille||2010| |Lamareille| 
et al. |2QQ4| ). Specifically, this diagram compares the ra¬ 
tios of equivalent widths [O III] A 5007/H/3 to [O II]/H /3 
and we can measure these ratios directly from the SOAR 
spectra. 

3.4. Reddening Correction 

HST and SOAR observations (when stated explicitly) 
are correc ted for intrinsic dust re ddening in a manner 
similar to McDonald et al. (2013|). Assuming continu¬ 
ous star formation, the young-stellar continuum L v (v) 
in the rest wavelengths 1500 — 2800A is flat, and we can 
thus construct reddening maps for the CLASH BCGs by 
assuming the gradual slope in flux we observe between 
pairs of UV images is due to mild dust extinction. In 
most cases we select the F275W and F336W filter to 
calculate the reddening. Filter selections for individual 
BCGs are given in Table 2. Even though this procedure 
assumes continuous star formation, the effect of SFH is 
relatively small- a starburst with a finite burst duration 
above 10 Myr will have a slope that differs from a model 
with continuous star formation by < 5%. 

We construct spatially resolved reddening estimates by 
covering the UV emission features in each BCG with grid 


squares 0.195” on a side, corresponding to 3x3 pixels 
in the drizzled images. Grid squares of this size are a 
compromise between sensitivity and spatial resolution, 
and make our corrections direct ly comparable to the red - 
dening correction employed in McDonald et al. (2013). 
Squares which do not have a minimum of 5 a of U V flux in 
the two filters used to calculate E (B — V ) are rejected, 
leaving behind a grid covering just the significant UV 
emission, like the example given in Figure [4] 

We ca lculate two redd e ning m aps for each BCG, one 
using a Calzetti et al. (2000) extinction curve, and 
one using t he Mi lky Way dust curve parameterized in 
O’Donnell (1994). When using the Milky Way dust 


model, we avoid the effect of the 2175 A bump by both 
subtracting the best-fit model of the bump in |Fitzpatrick| 


& Massa (1986) and by choosing UV filters that minimize 
coverage of the bump. In each grid square, E (B — V) is 
calculated by solving for the dust extinction necessary 
to flatten the slope between the UV-upturn-corrected 
fluxes in the two filters. The resulting reddening maps 
are Gaussian smoothed with a 0.195” kernel, in order to 
blur out the effects of binning. 

Out of the 20 clusters examined from the CLASH sam¬ 
ple, 11 have sufficient UV flux to estimate intrinsic red¬ 
dening over at least part of the BCG. Images of the esti¬ 
mated intrinsic reddening maps are shown in the right¬ 
most panel for each of these clusters in Figure [3] We 
find that, with few exceptions, E (B — V) < 0.5, which is 
consistent with the dust content typical of clu ster BCGs 
(|Crawford et al.||1999| |McDonald et al.||2011 ). 

Won-resolved dust extinctions for UV-faint BCGs and 
faint regions in UV-bright BCGs are calculated using 
F140W-IRAC colors assuming an underlying early-type 
stellar population. These bands are suitable for estimat¬ 
ing extinction in these galaxies and regions since the 
F140W band is slightly extinguished at the rest wave¬ 
lengths of CLASH BCGs while the IRAC bands are es¬ 
sentially reddening free. Furthermore, the F140W-IRAC 
color at these redshifts is insensitive to the SFH, so the 
assumed underlying population does not affect the re¬ 
sulting dust estimate. Spitzer/IRAC 3.6/im and 4.5/im 
fluxes are available for gill CLASH clusters except Abell 
1423 (Moustakas et al.||in prep.|). We use fluxes mea¬ 
sured in 3.0” diam eter apertures and apply the aperture 
corrections used in Sanders et al. (2007), selecting either 
the 3.6 /im and 4.5 /im band for each BCG separately in 
order avoid the polycyclic aromatic hydrocarbon feature 
at 3.3 /im rest-frame. 

Detections of H 2 vibrational modes are also pre valen t 


in the IR betwee n 5-25/im in starforming BCGs (Don- 
ahue et al.| 20ll|). We note that the presence of these 
Tines in IRAC Liters may affect our estimate of the dust 
reddening; however, we are mostly relying on this esti¬ 
mate of the reddening in BCGs with little evidence of 
ongoing star formation where we would not expect there 
to be wide vibrational H 2 lines. 

We use the spatially resolved extinctions to correct 
observed fluxes in UV-bright structures and the non- 
resolved extinctions to correct fluxes outside these struc¬ 
tures and in UV-faint BCGs. To correct line luminosi¬ 
ties, we adopt the re latio n E(B — V )* = 0.44 E(B — V) gas 
reported in Calzetti et al. (2000). While this is the em¬ 
pirically observed relation between the extinction of neb- 
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Figure 3. The reddening corrected UV luminosity map for CLASH BCGs are shown in the left panel. The region used to 
extract Luv for each BCG is outlined in black. Likewise, the Ha images are shown in the center panel. The estimated Calzetti 
reddening map is in the rightmost panel, in units of E(B — V). The uniform color code in the E(B — V ) ranges from 0.0 
(black) to 0.8 (red), with typical values between 0.1 and 0.5. White regions in the E(B — V) maps have insufficient UV flux to 
perform the reddening estimate. Galaxies shown are Abell 383, MACS0329, MACS0429, MACS1115, MACS1423, MACS1720, 
MACS1931, MS2137, RXJ1347, RXJ1532, and RXJ2129, in that order. 


ular and stellar emission in starburst galaxies, our choice 
of extinction model will introduce a systematic uncer¬ 
tainty, since starburst BCGs may differ from the star- 
burst galaxies used to calibrate this relationship. 

The reddening correction multiplies our values for L uv 



dening corrected UV SFR for RXJ1532 we find is 97 db 4 
M e yr -1 (only accounting for statistical uncertainty). 
For this BCG, the reddening corrected UV SFR yields a 


rate comparable to the IR estimate. 

3.5. Broadband Aperture Photometry 

To extract photometry for estimating the mean UV 
luminosities (L v ,uv) in UV-bright BCGs, we chose re¬ 
gions that contain all of the substantial UV flux with the 
aid of the ds^] contour tool (see Figure |3|. These re¬ 
gions trace surface brightness contours of 7.14xl0 24 erg 
s -1 Hz -1 pix -2 , corresponding to a star formation rate 
surface density of 10 -3 M© yr -1 pix -2 . For the ten UV- 
faint BCGs (those lacking sufficient UV flux to make a 
spatially resolved estimate of reddening), we measured 

2 http://ds9.si.edu 
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Figure 3. Continued 
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Table 2 

UV and H a Filters for BCG Luminosities 
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Cluster 

z a 

UV filters 

E(B-V) filters 

Ha ‘continuum’ 

Ha ‘line’ 

Abell 209 

0.209 

F225W, F275W 
F336W 

— 

F850LP 

F775W 

Abell 383 

0.187 

F225W, F275W 

F225w, F275w 

F850LP 

F775W 

Abell 611 

0.288 

F225W, F275W 
F336W 

— 

F775W 

F850LP 

Abell 1423 

0.213 

F225W, F275W 
F336W 

— 

F850LP 

F775W 

Abell 2261 

0.224 

F225W, F275W 
F336W 

— 

F625W 

F775W 

MACS0329.7—0211 

0.450 

F225W, F275W 
F336W, F390W 

F275W, F390W 

F775W 

F850LP 

MACS0429.6—0253 

0.399 

F225W, F275W 
F336W 

F275W, F336W 

F775W 

F850LP 

MACS0744.9+3927 

0.686 

F275W, F336W 
F390W, F435W 

— 

F140W 

F105W 

MACS1115.9+0129 

0.352 

F225W, F275W 
F336W 

F275W, F336W 

F775W 

F850LP 

MACS1206.2-0847 

0.440 

F225W, F275W 
F336W, F390W 
F225W, F275W 

F275W, F390W 

F775W 

F850LP 

MACS1423.8+2404 

0.545 

F336W, F390W 

F275W, F390W 

F125W 

F105W 



F435W 




MACS1720.3+3536 

0.391 

F225W, F275W 
F336W 

— 

F775W 

F850LP 

MACS1931.8-2635 

0.352 

F225W, F275W 
F336W 

F275W, F336W 

F775W 

F850LP 

MS2137—2353 

0.313 

F225W, F275W 
F336W 

F275W, F336W 

F775W 

F850LP 

RXJ1347.5—1145 

0.451 

F225W, F275W 
F336W, F390W 

F275W, F390W 

F775W 

F850LP 

RXJ 1532.9+3021 

0.363 6 

F225W, F275W 
F336W 

F275W, F336W 

F775W 

F850LP 

RXJ2129.7+0005 

0.235 

F225W, F275W 
F336W 

F275W, F336W 

F850LP 

F775W 

RXJ2248.7-4431 

0.348 

F225W, F275W 
F336W 

— 

F775W 

F850LP 

CLJ1226.9+3332 

0.890 

F336W, F390W 
F435W, F475W 

— 

F105W 

F125W 

MACS1311-0310 

0.499 

F225W, F275W 
F336W, F390W 

— 

F125W 

F105W 


a Redshifts are the sam e as those quoted in |Postman et ab| ( |2012| ), unless otherwise stated. 
b ( Crawford et al.|l999| > 



Figure 4. Illustration of the method for extracting reddening 
information from UV filters, for RXJ1532.9+3021. The final 
reddening map is shown in the right. The layout of grid 
squares used to calculated E(B-V) is shown overlaid on the 
F275W filter on the left. 


fluxes inside circular apertures that were as large as pos¬ 
sible but still excluded satellite galaxies. We opted for 
this strategy for measuring photometry because it makes 
it straightforward to capture the flux in the highly irreg¬ 
ular morphologies present in the CLASH BCGs. We cal¬ 
culated H<a+[N II] luminosities using the same apertures 
selected for significant UV flux. 

Line luminosities measured in our spectra were mea¬ 
sured in rectangular apertures that, while similar in area 
in most cases to the apertures we used to measure broad¬ 
band luminosities, nonetheless do not cover the same 
regions of the BCGs that were included in our calcula¬ 
tions of L uv and L hcx+[nii]- Therefore, when comparing 
line luminosities obtained from SOAR spectroscopy with 
broadband luminosities, it was necessary to measure the 
broadband luminosities in the apertures corresponding 
to the position and shape of the slit. Since our reduced 
spectra are 1-D, we used the mean value of E(B — V ) 
in each rectangular aperture to estimate the extinction 
correction for both the spectral line luminosities and the 
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UV and H<a+[N II] luminosities we compare them to. 

3.6. Stellar Population Properties 

We use iSEDfit to calculate the probability dis¬ 
tributions for model approximations of the SED in 
RXJ1532.9+3021, either in single apertures or in in divid ¬ 
ual pixels t o crea te stellar parameter maps. See [M ous- 
takas et al.|(|2013|) for details on iSEDfit. SEDs are com¬ 
posed of tiuxes extracted from identical apertures (or in¬ 
dividual pixels), in each of the 16 bands of CLASH HST 
photometry. We do not correct for fluxes for intrinsic 
reddening, since iSEDfit allows for local extinction to 
be treated as a fit parameter. 

When creating parameter maps using SEDs fit on in¬ 
dividual pixels, we PSF-match each image to the F160W 
PSF and extract SED for each pixel using the PSF- 
matched photometry. For approximately normally dis¬ 
tributed parameters, we assigned the mean values for the 
posterior probability density functions (PDF) for model 
SED parameters for each pixel SED to the locations of 
pixels in order to create maps of the BCG. The model 
predicted rest-frame flux for H<a+[N II] provides a san¬ 
ity check on the physical parameter maps for the 16 
band SED, since we can compare its morphology to the 
L Ha+[Nii] ma P described in Section 3.2. The two im¬ 
ages are shown in Figure [5] with matching coordinates, 
where it can be seen that the filamentary features in the 
L hol+[nii\ ma P correspond to the features derived from 
the SED. The two images depict the same pair of Ha 
‘bulges’ in the center of the BCG as well. 





Figure 5. Left: Map of the Ho+[N II] luminosity constructed 
using iSEDfit. Right: Ha+[N II] luminosity map estimated 
by scaling F850LP and F775W images. We juxtapose the two 
to demonstrate the morphological similarity b etw een them. 
The structures in our starburst maps in Figure [I5| appear to 
be dominated by the filamentary structure in the Ha+[N II] 
image. 


4. RESULTS 

4.1. Broadband Luminosities and Star Formation Rates 

Mean UV and H<a+[N II] luminosities are given in Ta¬ 
ble 3. We present the reddening corrected ( Luv) both for 
a Calzetti reddening law and for Milky Way-type dust. 
Luminosities are converted to SFRs using the Kennicutt 
SFR calibrations. There are several sources of poten¬ 
tial scatter in this estimate, including contamination of 
the UV luminosity by AGN activity, and va riatio ns in 
the IM F and SFH of the stellar population (Kennicutt 
(1998) assumes a continuous SFH and a Salpeter ( |1955p 



Figure 6. L uv and L hol+[nii] are plotted for the regions 
tracing UV emission. The luminosities are corrected for dust 
obeying a Calzetti reddening curve. The light blue band 
shows the region corresponding to where the two luminosi¬ 
ties predict the same SFR according to the Kennicutt SFR 
calibrations, to within the 0.3 dex scatter in the Kennicutt 
UV calibration. Absent error bars are too small to see as 
plotted. 


IMF). For the L uv based SFR estimates, we have not 
attempted to correct for these effects. However, the UV 
features we observe are not likely to be due to AGN ac¬ 
tivity, which is ruled out by the co mplic ated UV mor¬ 
phology reported in Donahue et al.| (|2015|). _ 

In o rder to estimate Ha based SFRs using [Kennicutt | 
( 1998| ), we need to estimate the ratio [N Il]/H<a. The line 


ratio of [N II] to Ha can vary between BCGs, and within 
filamentary structures in BCGs, so whatever choice we 
adop t will add scatter to our estimate of the Ha based 


SFR (McDonald & Ve illeux|2QQ9 McDonald et al. 2014a 


Crawford et al. 1999). ~This ratio is typically 0.5 for opti 
cal galaxies (|Kennicutt|l992||Kewley et al. (2001b, 2004) 
However, the ratio |JN II]/Ha is often larger thah~TT 5 
in BCGs (Heckman et al. 1989b Crawford et aLpf)^ 


For Ha luminous BCG s in the X-r ay selected samp 


of clusters in jCrawford et al. (1999), the typical [N II] 
A6584/Ha is 1.1 ± 0.4. We adopt this value in order to 
calculate SFRs using L Ha+[Nii]i bearing in mind that 
this is a rough approximation with considerable scatter. 
Nonetheless, we believe that variation in [N II] /Ha is a 
secondary consideration for the purp oses o f estim ating 
SFRs, given the scatters in the Kennicutt| (fT998) cali¬ 
brations between L# a and SFR and between Ljjy and 
SFR. 

The correlation between UV and H<a+[N II] luminosi¬ 
ties is shown in Figure |6l The two luminosities are 


broadl y cons istent with the ratio expected from Ken- 
~ ([1998}, in that the SFR estimates derived from 


nicutt 


the UV and H<a+[N II] luminosities are consistent with 
each other. This result is d iffers from the findings of Mc-| 
Donald et al.| ( |2Q10| 2011), since they find on average 
the UV/Ha ratio is slightly lower than that predicted 
from the Kennicutt relationships. This is most likely be- 



















































































Table 3 

UV and H a Filters for BCG Luminosities 
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Cluster 

10 27 

L uv 
(Cal) a 
ergs s —1 Hz -1 

L[/V 

(MW) 6 

^Hoi+[NII] 

(Cal) a 
10 41 ergs s -1 

UV SFR“ 

M© yr” 1 

Ha SFR a 

Area c 

kpc 2 

Abell 209 

0.1 ±0.3 

0.1 ±0.3 

- 

0.01 ±0.04 

- 

- 

Abell 383 

21.4 + 2.3 

37.2 + 3.3 

5.9 + 1.0 

3.0 ±0.3 

1.9 ±0.6 

14.26±0.02 

Abell 611 

< 0.3 d 

< 0.3 

- 

< 0.04 

- 

- 

Abell 1423 

0.5 ±0.2 

0.6 ±0.3 

- 

0.07 + 0.03 

- 

- 

Abell 2261 

0.2 ±0.2 

0.2 ±0.2 

- 

0.02 + 0.02 

- 

- 

MACS0329.7—0211 

302.3 + 12.1 

349.9 + 12.6 

254.0 + 24.1 

42 + 2 

80 + 21 

173.9±0.2 

MACS0429.6—0253 

200.1 ± 11.2 

234.2 ± 10.6 

103.3 + 11.8 

28 ±2 

33 + 9 

72.2 ±0.2 

MACS0744.9+3927 

4.2 ±0.6 

4.6 ±0.7 

- 

0.6 ±0.1 

- 

- 

MACS 1115.9+0129 

91.6 + 6.5 

92.8 + 5.5 

10.7+1.7 

13 ± 1 

3.4 + 1.0 

52.32 ±0.05 

MACS 1206.2—0847 

20.8 + 4.4 

25.4 + 3.9 

- 

2.9 ±0.6 

- 

- 

MACS 1423.8+2404 

193.7 + 8.2 

212.4 + 8.3 

103.3 ± 11.1 

27+1 

33 + 9 

121.2 + 0.2 

MACS1720.3+3536 

7.6 ±1.5 

22.5 + 5.7 

13.8 ± 1.6 

1.1 ±0.2 

4.4 ± 1.2 

10.08 + 0.02 

MACS1931.8—2635 

1975 ± 135 

1756 ± 79 

422.1 ±49.5 

280 ± 20 

130 ± 40 

331.5 ±0.5 

MS2137-2353 

43.7 ±6.7 

22.4 + 2.2 

35.3 ±9.9 

6.1 ±0.9 

11 ±4 

15.41 ±0.07 

RXJ1347.5—1145 

153.7 + 9.0 

175.4 + 8.6 

36.7 + 6.9 

22 + 1 

12 ±4 

71.3 + 0.2 

RXJ1532.9+3021 

694.0 + 25.0 

765.6 ± 24.5 

429.9 + 29.1 

97 + 4 

140 ± 40 

308.1 + 0.2 

RXJ2129.7+0005 

12.4 ±6.9 

9.3 ±3.0 

32.5 ±31.1 

1.7 + 1.0 

17.2 ± 16.4 

10 ± 10 

RXJ2248.7—4431 

0.8 ±0.4 

0.9 ±0.4 

- 

0.1 ±0.1 

- 

- 

CLJ1226.9+3332 

6.2 ±0.4 

6.5 ±0.5 

- 

0.9 ±0.1 

- 

- 

MACS1311.0—0310 

3.9 ±2.9 

4.0 ±3.0 

- 

0.5 ±0.4 

- 

- 


a Calzetti model dust was used to calculate the reddening correction. 
b Milky Way model dust was used to calculate the reddening correction. 

c Areas of UV emitting regions as observed with the F336W filter. Uncertainties are calculated by using Monte Carlo draws to sample the 
distribution of the number of pixels containing positive flux in the regions shown in Figure [3] 
d 3cr upper limits are shown. 


cause they do not correct for extinction due to dust in 
the BCG. Indeed, they propose adding a correction of 
E (B — V) = 0.2 to their data, which would make their 
results consistent with continuous star formation, and 
this value is typical for the dust extinction we observe in 
CLASH BCGs. 

UV SFRs are correlated with the areas of the star 
forming region in CLASH BCGs (Figure [7]). CLASH 
BCGs have an average SFR surface density ((ESFR)) of 
~ 0.3 M 0 yr -1 kpc -2 , with typical values ranging be¬ 
tween ^ 0.1 — 0.4 M 0 yr -1 kpc -2 . Areas of the UV flux 
emitting regions were measured using the F336W filter. 
In order to calculate the uncertainty, we sampled the dis¬ 
tribution of fluxes in each pixel in the regions shown in 
Figure [3] using a Monte Carlo method, which we used 
to create a distribution of flux-emitting areas. The ex¬ 
ception to this is MACS1931.8—2635, which exhibits a 
(ESFR) of 0.83 =b 0.06 M 0 yr -1 kpc -2 . 


4.2. SOAR Spectra Results 

The CLASH BCG UV and [O II] luminosities, dis¬ 
played in Figure |8j scale with each other, but produce di¬ 
vergent SFR estimates. However, UV and R/3 luminosi¬ 
ties (Figure [9]) have a tight correspondence and produce 
consistent estimates of the SFR in CLASH BCGs. The 
agreement between UV and H/3 based SFRs is tighter 
than the agreement between UV and H<a+[N II] derived 
SFRs, which is to be expected considering the limited 
precision of SFRs estimated using broadband H<a+ [N II] 


The SFR-L[o//j relation we use is calculated in Kew- 


ley et al. (2004) by using a sample mean [O Il]/H<a to 



Figure 7. The UV derived SFR is shown as a function of the 
area of the UV flux emitting regions measured in the F336W 
filter. The dotted line represents an (ESFR) of ~ 0.3 M 0 
yr -1 kpc -2 . Absent error bars are too small to be seen as 
plotted. 


convert from the Kennicutt (1998) SFR-L# a relation to 
an SFR-L[o/j] relation. However, the theoretical value 
of L[o//]/SFR depends on the metallicity of the nebular 


region (peaking near Z ~ 0.5 Z 0 ) as well as ionization 
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Figure 8. The distribution of BCG UV luminosities, L uv, 
compared to [0 II] luminosities, L[on], for all BCGs with 
SOAR coverage of [0 II]. UV luminosities were measured 
in rectangular apertures that correspond to the spectral slit 
placements. The light blue band corresponds to where these 
luminosities predict the same star formation rate to within 
0.3 dex, as in Figure [6] 


para meter (peaking near q ~ lx 10 7 cm s x ) (Kewley 
et al.| 2004). CLASH BCG redenning corrected flux ratios 


[O III] / [O II] are typically ^0.1, implying an ionization 
parameter ne ar q ~ 1 x 10 7 cm s" 1 for solar and sub-solar 
metallicities (jKewley & Dopita||2002|). The combination 
of these two parameter depedencies may explain the sys¬ 
tematic tension between UV and [O II] SFRs. Further¬ 
more, the offset between UV and [O II] SFRs we observe 
in CLASH clusters is consistent with the observation in 


Kennicutt (1998) that L[o/q/SFR is typically boosted 
in sfarbursts relative to galaxies undergoing continuous 
star formation by a factor of > 2. 

SOAR spectra were also used to constrain the source 
of the photoionizing emission we observe. We place ac¬ 
tive CLASH BCGs on the blue-line diagnostic diagram 
for distin guishing star bur sting g alaxi es from AGN, de - 
scribed in Lamareille et al. (2004) and|Lamareille (2010), 
as well as on the HFT' diagram, in f igure TO We 
cannot directly separate Ha from [N II] in our broad¬ 
band Hct+[N II] fluxes, so when available, we use line 
fluxes from the SDSS Data Release to determine 
the locations of C LASH BCGs on the BPT diagram 
(Alam et al.||2015). In order to place the remaining 
active BCGs on the BPT diagram, instead of compar¬ 
ing [O III] /HyS to [N Il]/Hcq we compare [O III] /H/3 

to X = 0.75 (~ l)- Our expression for X as¬ 
sumes case B recombination. Regardless of the presence 
of AGN emission, case B recombination allows us to de¬ 
rive a reasonable estimate of the ratio of Ha to H/3, since 


for svstems with hvdrogen densitie 

is in the range 10 3 —10 6 

cm 3 , Hct/H/3 ~ 2.7-3.2 ( 

Netzer 

2013). Because AGN 


often produce harder photoionizing spectra than young 


Figure 9. BCG UV and H [3 luminosities for all BCGs with 
SOAR coverage of H/3. UV luminosities were measured in 
rectangular apertures that correspond to the spectral slit 
placements. In order to determine the region in the plot where 
the two luminosities predict consistent continuous SFRs, we 
scaled the Kennicutt law relating L Ha to the SFR by a factor 
of 2.85. The light blue band is analogous to thes band de¬ 
picted in Figures ISland [8] Typical uncertainties on L up are 
~10%, and typicaTuncertainties on L uv are ~ 5%. 


stellar populations, our assumption will tend to bias X 
slighter higher than [N Il]/H<a in the presence of an AGN. 
However, such a bias will cause the estimated line ra¬ 
tios to appear more ‘AGN-like,’ so the resulting BPT 
diagram is a conservative estimate of the contribution 
of ongoing star formation to the line ratios we observe. 
The positions of BCGs on the BPT diagram determined 
using X depend on both the accuracy of our reddening 
corrections and broadband Hct+[N II] estimates. There¬ 
fore, their value is primarily as a consistency check of the 
blue-line diagram. 

For the blue-line diagram, we use equivalent widths ob¬ 
served with SOAR for all CLASH BCGs except Abell 383 
and MACS1423.8+2404. Since these BCGs have incom¬ 
plete line flux data from the SOAR spectra (Abell 383 
does not have an [O II] measurement and MACS1423 
has an upper limit for [O III] A5007 estimated from an 
upper limit of [O III] A4959) but were observed in Data 
Release 12, we use SDSS equivalent widths for these lines 
instead. In general, our SOAR spectra are better suited 
to observing extended nebular emission in CLASH BCGs 
because we were able to place the slit to maximize cover- 


galSpecn galaxy sample (Brinchmann et al. 

to 

o 

o 

Kauff- 

mann et al. 2003a Tremonti et al. 2004). C 

rASITBCGs 


low [O III]/H/5 and high [O II]/H/3 relative to the SDSS 
dataset. Our results imply that most of the BCGs lie 
i n the co mp osite s tar forming-LINER region described in 


Lamareille (2010), with the exception of MACS1931.8- 
KT" 


3 http://drl2.sdss3.org/ 


4 http://www.sdss.org/drl2/spectro/galaxy_mpajhu/ 
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Figure 10. Left: The blue line diagnostic diagram for line-emitting CLASH BCGs. Individual BCGs appear labelled on the 
plot. Galaxies from galSpec are shown color-coded by BPT classification. Green galaxies are starforming galaxies, light blue 
galaxies are composites, blue galaxies are LINERs and orange galaxies are AGNs. The positions of MACS1931 and RXJ1532 
are highlighted, since these two are the largest starbursts in the CLASH BCG sample. Black points are measured from SOAR 
data, while red point are measured from SDSS data. Right: The BPT diagram incorporating CLASH broadband line flux 
estimates. Red points are measured from SDSS data, while black points are estimated using a combination of SOAR fluxes and 
HST broadband Ha+[N II] flux estimates. For the black points, we defined a proxy for [N II]/Ha on the x-axis of the BPT 

diagram, A, to be 0.75 ( — l)- K e gi° n labels match color coding and approximate locations of regions populated by 


1.8577/3 

different classifications of galaxy. 


The BPT diag ram shows the Kewley et al. (2001a) 
line in blue and Kauffmann et al. (2003b) line in green. 
CLASH BCGs are distributed in the starforming and 
composite-starforming regions of the diagram. The 
CLASH BCGs cluster around log([0 IIlJ/H/3) ~ —0.3, 
which puts them below the BPT discriminating bound- 
ary between starforming galaxies and AGN (| Kauffmann 
et al. 2003b). The exception to this is MAC Si 931, which & 


is consistent with emission powered predominantly by 
star formation. We observe an X-ray AGN in the Chan¬ 
dra image of MAC SI 931; however, given the extent of the 
UV emission region it makes sense to classify the BCG 
as starbursting. For the most part, the blue-line and 
BPT indicate consistent sources powering line emission 
in CLASH BCGs. 

Based on these diagnostics, we conclude that the line 
emission in most of the BCGs is either predominantly due 
to ongoing star formation, or to a composite starforming- 
LINER-like source. In particular, the two most UV lu¬ 
minous BCG in our sample, MACS1931 and RXJ1532, 
are consistent with star formation being the main pho¬ 
toionization mechanism when taking into account both 
diagnostic diagrams. MACS1720, MS2137 and possibly 
Abell 383 may be LINERS, although much of their UV 
and H<a+[N II] flux is not in a nuclear emission region. 
Likewise, the majority of the composite galaxies fall into 
the composite starforming-LINER classification. This is 
consistent with previous resu lts finding LINER-like_emis- 
sion in cool-core B CGs ( Veron-Cetty fe Veron|200Q Ed- 


wards et al.||2007 ). However, while the emission line di¬ 


agnostics in these BCGs are LINER-like, they cannot be 
LINE Rs since they cannot be powered by a central black 
hole (Heckman et al. 1989b). Several hypotheses have 


been proposed for the source of this extended LINER- 
like emission. Stellar populations may be responsible 
for this emission, which could be due to photoionization 


([Shields 

1992 

Olsson et al. 

2010 Loubser & Soechting 

20 Ufa b 

. Emission lines may also be due to nebular gas 


fc Voit|1991| Werner et al.) 2013), 

(Begelman & Fabian 1990), or co 

et al.||l989[ [herland et al.|[2009 ). 


, turbulent mixin g layer s 


collisional heating (Sparks 


4.3. Correlation with ICM X-Ray Properties 
4.3.1. Core Entropies 

SFRs derived from Luv in our sample are correlated 
with the X-ray properties of the ICM. Figure [TT] shows 
the relationship between CLASH BCG SFRs and ICM 
core entropies. The core entropy Kq used in the present 
study is defined to be the innermost bin of the entropy 
profile in ACCEPT, and is a proxy for t he existen ce of 
a cool core in a galaxy cluster (Hudson et al.|2010). En¬ 
tropy as measured by X-rays is defined to be 


K = kT x n ~ 2 / 3 , 


( 2 ) 


where T x is the X-ray temperature in keV, n e electron 
density in cm -3 , and k is the Boltzmann constant. 

Low values of Kq typically accompany activity in 
BCGs. BCG activity, such as elevated NUV flux rela¬ 
tive to the predicted quiescent UV emission, is observed 
to occur only in clusters where Kp is < 30 KeV cm 2 
dCavagnolo et al." 2008 Hoffer et al. 2012). McDonald 
et al.| (|2010|) reported on this phenomenon as well with 
resolvea Ha/ emission maps in low-redshift BCGs. In 
Donahue et al. (2015), a similar entropy threshold was 
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found for the UV-NIR color of CLASH BCGs, indicating 
the threshold does not change substantially out to z^0.5. 

Here, we demonstrate that a tight correlation exists 
between reddening corrected SFRs and Ko- Specifically, 
all of the BCGs with an SFR > 10 M 0 yr -1 have a core 
entropy consistent with a value < 30 KeV cm 2 (see Fig¬ 
ure UU. Meanwhile, BCGs that lack significant UV or 
HaH-[JNI II] luminosities occupy a range of core entropies 
that extends up to ^ 200 keV cm 2 . Considering the 
difference between the two observables, the tight corre¬ 
spondence between SFR and core entropy is compelling. 



Figure 11. BCG SFRs compared to cluster core entropies. 
Core entropies below ~ 30 KeV cm 2 are associated with the 
onset of UV emission in BCGs, and all of the CLASH BCGs 
with a UV SFR > 10 M 0 yr -1 reside clusters with a core 
entropy < 30 KeV cm 2 . The grey vertical line depicts this 30 
KeV cm 2 entropy threshold. Core entropies are taken from 
the ACCEPT cluster profile archive. For several points, er- 
rorbars for the SFR and Ko are too small to be depicted on 
this plot. 


4.3.2. Core M g Estimates 


ICM cooling rates inferred by Chandra observations 
of clusters bear little relation to the actual cooli ng rat e 
in cool-core clusters (e.g McNamara & 0’ConneII|[l989j). 
However, we may to be able to find a relationship Be¬ 
tween star formation and a simple proxy for cooling at 
radii where we hypothesize ICM cooling actually occurs. 


The quantity we examine, M^ (r), is analogous to the pre¬ 
dicted cooling rate as a function of radius, and is defined 
by 


M, (r) 


47r So Pg (?) r 2 dr 

tcool (j *) 

M enc (r) 

tcool (0 


(3) 


where p g (r) is the azimuthally average X-ray gas density, 
and t coo i (r) is the averaged X-ray cooling time at radius 


r. M enc (r) is the gas mass enclosed in the radius r. We 
use 

+ 3 nkT 

tco ° l = 2n e n„A (Z,T) 

to define the cooling t ime, and use the assumption in 
Cavagnolo et al. (2009) that n ~ 2.3 nn • The cooling 
c urve A (Z,T) i s estim ated by interpolating the Suther¬ 
land & Dopita (1993) cooling function at solar metallic- 
ity. 

We choose to measure M^ (r) at r = 35 kpc (M^^s), 
and at the radius in each cluster where the ratio be¬ 
tween the cooling time and free-fall time is t coo i/tff = 20 
(M^ 2 o). We calculate free-fall times by estimating clus¬ 
ter density pro files as the sum of N FW profiles derived 
from lensing in Merten et al. (2015) and singular isother¬ 
mal sphere protiles ot B CG stellar density derived from 
stellar mass estimates in Burke et al. (2015). The impact 


of the stellar mass component on our overall result is not 
substantial; however, we include it for completeness. 

M ^ r 35 is useful to measure because 35 kpc is typical 
of the maximum radius we observe H<a+[N II] and UV 
structures in CLASH clusters, and because this radius 
maximizes the correlation between UV SFR and M^ (r) 
(see Figure 12). This quantity can be calculated for all 
the CLASH clusters using ACCEPT data. 



100 150 200 250 

R [kpc] 


Figure 12. Pearson correlation coefficients are shown be¬ 
tween log(SFR) and log(M^) as a function of radius R. Green 
line denote the minimum correlation to rule out the null hy¬ 
pothesis at P<0.05. 


The choice of M^^o reflects the finding that BCG ac¬ 
tivity occurs in cl usters with a minimum t coo i/tff be¬ 
tween 4 and 20 (Voit et al. 2015). Since BCG activ¬ 
ity is associated with potentially cooling ICM gas where 
tcooi/tff < 20, M s . £20 measures the predicted cooling 
rate in gas that we suspect is directly involved in cooling. 
Radii where t coo i/tff = 20 are listed for each cluster in 
Table 4. Clusters where t coo i/tff > 20 in the innermost 
bin of the cooling time profile calculated from ACCEPT 
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data are not included. 


Table 4 

tcooi/tff Threshold Radii 


Cluster 

Radius [kpc] 

Abell 383 

50 + 5 

MACS0329.7—0211 

79 + 5 

MACS0429.6—0253 

83 + 9 

MACS0744.9+3927 

50 + 11 

MACS1115.9+0129 

108 ± 17 

MACS1423.8+2404 

83 + 8 

MACS1720.3+3536 

68 + 6 

MACS1931.8-2635 

105 ± 14 

MS2137—2353 

83 + 6 

RXJ 1347.5—1145 

69 + 18 

RXJ 1532.9+3021 

121 ± 12 

RXJ2129.7+0005 

59 + 8 

MACS1311.0-0310 

83 + 24 


We show the relationship between the UV SFR and 
both M ^ r 35 and M^o in Figure 13| The dashed lines 


denote, from right to left, where rue BCG is forming 
stars at 100 %, 10 %, 1 %, and 0 . 1 % of the cooling rate 
implied by M^. If we interpret the SFR as a proxy for the 
actual cooling rate in this system, and interpret as the 
‘potential’ cooling rate in the absence of feedback, then 
the larger starbursts are cooling much more efficiently 
than smaller starbursts, and these lines indicate where 
the ‘efficiency’ is 100 %, 10 %, 1 %, and 0 . 1 %. 

We fit trend lines to the data for both the SFR-M ^ r 35 


and SER-M^^o relations using orthogonal least squares 
regression. We find that the slope on the trends fit to the 
two datasets (0.35+0.05 for SFR-M P;r 35 and 0.27+0.06 
SFR-M^o) are nearly consistent, leading to the con¬ 


clusion that the two definitions of measure a similar 
quantity. 

Several limitations impacting our measurements may 
affect how tightly correlated SFR-M^ (both M^r -35 and 
M^o) appear to be in our data. Gas density profiles 
in ACCEPT have a limited resolution (between 10—30 
kpc per bin depending on the CLASH cluster), so values 
for M enc are typically calculated by interpolating on the 
central few bins of each profile. Temperature profiles are 
less well resolved than p g profiles, which adds scatter to 
our estimate of t coo i (r) profiles. Deeper X-ray observa¬ 
tions will beat down the systematics in M^, and a larger 
sample of cool core clusters will allow us to more pre¬ 
cisely constrain the SFR-M^ relationship and examine 
the effects of sample selection. With our current data, 
we establish that a relationship exists between these two 
quantities, and that this relationship implies that as the 
BCG SFR increases, there is a steady increase in the ra¬ 
tio of ongoing star formation relative to the predicted 
cooling time in the reservoir of hot gas. 


4.4. Model Fitting to RXJ1532.9+3021 

We adopt RXJ1532.9+3021 as a case study and use 
its HST SED to delve into the SFH of the BCG in this 
galaxy cluster. RXJ1532 exhibits the second highest star 
formation in our sample and is replete with UV and Ra 





Figure 13. Relationship between SFR and measured at 
35 kpc (top) and radii where t coo z/t// =s 20 (bottom). Solid 
lines fit the data. Dashed lines indicate (from furthest right) 
where SFR is 100% of M^, 10%, 1%, and 0.1%. 


bright filaments and knots. This BCG makes a better 
case study than the strongest star forming galaxy in our 
sample, MACS1931.8-2635, because the latter exhibits 
a strong X-ray AGN which could complicate pixel-scale 
SED fitting. RXJ1532 also has detailed auxiliary data 
including an SDSS spectrum c overing Ra an d [N III and 


a deep Chandra observation (Hlavacek-Larrondo et al. 


2013). RXJ1532 is therefore an excellent prototype for 


exploring the characteristics of the star forming regions 
in CLASH BCGs. 


We both fit a single-aperture SED from the region 
used to calculate L uv and constructed stellar parame¬ 
ter maps. The best-fit SFR in the single-aperture SED 
is II 8 I 42 5 M© yr -1 . The average extinction is E(B-V) 
= 0.27 ± 0.07. Since the posterior probability distribu- 
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Figure 14. We show the two dimensional posterior probability distributions for A t^xSFR, A Burst Mass, A t^x Total Mass, 
SFRx Burst Mass, SFRx Total Mass, and Burst Massx Total Mass. Black points denote the individual Monte Carlo draws used 
to construct the posterior probability distribution, and contours are the 68%-95%-99.7% contours calculated from the kernel 
density estimate smoothed distribution. The marginal distributions for individual parameters are shown by the histograms on 
the diagonal. 


tion of the SFR is close to log-normal, the best-fit value 
we report is the mode of the distribution and the un¬ 
certainty we report is the 68.3% confidence interval. The 
best-fit starburst is relatively long-lived, with a burst du¬ 
ration logA4 [Gyr] = —0.16 ±0.47, and massive, with a 
total burst mass log M b [M 0 ] = 11.03 ± 0.36. The burst 
parameters in the fit have degeneracies (see Figure 14), 
although the peaks in the probability distribution sug¬ 
gest that we are constraining the parameters. 

We show stellar parameter maps in Figure [l5j includ¬ 
ing A tb, the SFR, and the spatial distribution of the to¬ 
tal mass of the BCG. Star formation is concentrated in 
knots in the center of the BCG, along with a network of 
six bright filaments and several dimmer filaments. The 
SFR morphology is consistent with the Ha and UV mor¬ 
phology. The sum of the pixel SFR modes in the single¬ 
aperture SED region is 119 M© yr -1 , which matches the 
single-aperture value remarkably well. 

Two filaments point south (A in Figure p±| , two point 
northwest (B, Cl), one points southeast (C2), and one 
extends northward before bending to the east (D). The 
filament pointing southeast and the brighter nothwest 
pointing filament (Cl and C2) appear to lie along a single 
axis. The clumps in the core and the filaments in the 


north account for the bulk of the ongoing star formation. 
The starbursts in these structures are long-lived, on the 
order of 10 8 — 10 9 yr, which is consistent with the results 
of the single aperture SED. The filaments to the south, 
and the two that lie along a single southeast-northwest 
axis passing through the center of the BCG, are more 
than an order of magnitude younger than the rest of the 
burst system. 

The peak and western ‘bulge’ in the total mass cor¬ 
responds to the central knot morphology in the BCG. 
Meanwhile, two roughly conical ‘drop-offs’ appear to the 
east and west, and both of them extend outward from 
the positions of star-forming filaments. Similar drop-offs 
are be visible at the positions of the young, southern fil¬ 
aments. These drop-offs may reflect real deficits in the 
stellar surface density of the BCG, but may also be a 
consequence of the dust geometry in this system. Since 
dusty filaments in the BCG screen the elliptical stellar 
populations behind them, it is possible that the drop-offs 
are regions where the mass estimate is biased low due to 
the positions of filaments along the line of sight. 

In the Appendix, we demonstrate that the iSEDf it de¬ 
rived values for SFR and E(B — V ) agree with the broad¬ 
band values and show that the SED predicted H<a±[N II] 
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Figure 15. Two dimensional maps of SFR, log A t&, and the total mass surface density for RXJ1532. The bottom right panel 
provides labels for the morphological features in the parameter maps. The SFR map is in units of Mq yr _1 , log A t b in log Gyr, 
and the total mass in units of M© pix -2 . All values are expectation values, log A t b is masked to only depict regions where the 
likelihood of a starburst is > la. The pair of filaments to the south of the center of the BCG, and the filament along the NW/SE 
are noteworthy for being by > 1 order of magnitude younger than the average age of the starburst in this system. Possible 
conical shaped underdensities are roughly along the positions of bright filaments. At the redshift of the cluster, z=0.363, the 3” 
reference scale in the figures corresponds to about 15.2 kpc. 


line emission feature matches the SDSS spectrum. We 
also discuss the importance of the Ha feature to charac¬ 
terizing the SFH using SED fitting. The particulars of 
the SFH model (a uniform starburst superimposed on an 
exponentially decaying SFH) and our choice of parame¬ 
ter space are documented in the Appendix as well. 


5. DISCUSSION 

Half (10 out of 20) of the X-ray selected sample of 
CLASH clusters show evidence for significant (> 5 a) 
rates of reddening-corrected star formation using both 
UV and Ha indicators. CLASH BCGs occupy regions 
of line diagnostic diagrams that are typical of compos¬ 
ite starforming-LINER galaxies, and in several cases line 
emission may be primarily powered by star formation. 
This rate of incidence is substantially higher than pre¬ 
vious published rates of incidence of star formation or 
line emission in X-ray selected cl uster BCGs, which are 


in general closer to 2 0-30% (e.g. Crawford et al. 1999 
Edwards et al.||2007|). However, the CLASH X-ray se¬ 
lected sample differs from these populations of galaxy 
clusters, since it is comprised of high gas temperature 
(kT^ > 5 keV) clusters chosen according to a relaxation 
criterion based on X-ray morphology. The CLASH sam¬ 
ple of BCGs has an incidence of line emission similar 
to the incidence of line emission in REXCESS cool core 


clusters (70%) (Donahue et al. 2010). Our sample char¬ 
acteristics differ from REXCESS in that it is at higher 
redshift [z — 0.2 — 0.7 compared to z = 0.06 — 0.18) and 
along with being X-ray selected, CLASH clusters were 
selected for exhibiting relatively condensed, round X-ray 
isophotes. 

The trends between Luv and L Ha+[Nii] and between 
Luv an d L h /3 suggest applying the Kennicutt SFR cali¬ 
brations produces consistent star formation rates. How¬ 
ever, SFRs predicted using L [ 0 u] are systematically el¬ 
evated relative to UV based SFRs. This is in contrast 
to our findings using L Ha+[Nii] and L h/3, both of which 
predict SFRs consistent with Luv- These results may 
not be unusual for starburst galaxies, although they may 
also indicate that the [O II] emission line is being par¬ 
tially powered by an additional source heating the ionized 
gas. 

The SFRs in several BCGs are very large. In particu¬ 
lar, two galaxies exhibit SFRs > 100 M© yr -1 , and an 
additional five have SFRs > 10 M 0 yr -1 . The strongest 
star formers (MACS1931, RXJ1532) are forming stars 
several times more slowly than the Phoenix cluster, 
which to date exhibits the large st known SFR of a BCG 
dMcDonald eT!I1|lH3l |2014b[ ). However, MACS1931 
is noteworthy because its UV SFR (280 M 0 yr -1 ) is ~ 
40% of the cooling rate in the absence of heating (^700 
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M© yr -1 ) measured in Ehlert et al. (2011). The Phoenix 
SFR is ~30% of the cooling rate in the absence of heating 
measured in McDonald et al. (2014b), so it is plausible 
that MAC SI 931 and the Phoenix cluster harbor BCGs 
undergoing similar feedback events. The presence of an 
X-ray AGN in each BCG also suggests that the AGN is 
undergoing a similar evolutionary phase. Furthermore, 
MAC SI 931 is forming stars more densely than the rest 
of the CLASH sample, which suggests it is an outlier 
relative to other starforming BCGs. 


5.1. BCG-ICM Interactions 

Examination of core entropies implies that the ex¬ 
tended star forming features in CLASH BCGs are likely 
due to an interaction between the BCG and the envelop¬ 
ing ICM. Reddening corrected SFRs obey t he 30 keV 
cm 2 core entropy threshold reported in e.g. [Hoff er et al. 
(2012) - all the strong star formers (SFR > 10 M© yr _i ) 
fall at or below the threshold. From these results we con¬ 
clude that ongoing star formation in the BCGs is corre¬ 
lated with the thermodynamics of the surrounding ICM. 
It is plausible that a low ICM core entropy is necessary 
for the onset of star formation in these BCGs. However, 
it does not directly trigger star formation, as evidenced 
by the existence of low-SFR BCGs with core entropies 
below 30 keV cm 2 . 

We also analyze observables related to cooling in the 
low-entropy ICM surrounding BCGs, in order to better 
understand the interaction between the low-entropy ICM 
and BCG starbursts. We define two quantities, M ^ r 35 
and M g j 2 o, which approximate the cooling rate of ICM 
gas in the vicinity of the BCG. For both definitions of the 
cooling rate, we observe a similar trend between SFR and 


(2012), implementing bi-polar AGN jet feedback, fueled 
by cold gas, in an ENZO adaptive mesh simulation. They 
found that AGN feedback balanced cooling robustly, and 
generated episodic behavior on timesca les of a G yr, in 


M^. The positive correlations between SFR and both 

M^ r 35 and M. g j 20 are reasonable since the low-entropy 
gas near the BCG is a prime candidate for the reservoir 
of gas that cools to become star forming molecular gas. 
Since the correlation between SFR and M^ (r) drops as 
Mp (r) is measured at larger radii, these findings are con¬ 
sistent with the tension between observed ICM cooling 
rates those predicted by measuring ICM gas masses with 
t coo i less than a Hubble time. 

Recent theoretical work has made significant progress 
in identifying some of the key elements of the AGN feed¬ 
back processes in clusters of galaxies. Simulations in 
which cold g as drives the accret i on onto an AGN hav e 
been done by Gaspari et al. (2012); Li & Bryan (2014a|b); 
Li et al. (201oa|). Gaspari et al.|(2012|) modeled bi-polar 
AGN jets he ating an ICM atmosph ere, f ollowing up on 
the w ork of McCourt et al. (2012) and Sharma et al. 
(|2012|) , who found that circumgalactic gas in which heat¬ 
ing balances cooling gl obally becomes the rmally unstable 
when t coo i/iff < 10- Gaspari et al. (2012) found that this 
criterion also demarcated the transition to a multiphase 
medium in si mulations relying on AGN jet heating fueled 
by cold gas. Gaspari et al. (2013, 2015) extended this 
work by focusing on on the sub-pc scale behavior and 
demonstrated that chaotic cold accretion could signifi¬ 
cantly increase the accretion rate onto the central black 
hole, relative to the Bondi rate one would i nfer fr om the 
thermal state of the hot phase alone. Li fe Bryan (|2014a|) 
examined a mechanism similar to that in Gaspari et al. 


alignment with the work of Gaspari et al. (2012). The 
scale and filamentary nature of the cold filaments in Li 
& Bryan simulations resemble those we observe, and the 
ages of star forming filaments we recover in RXJ1532 are 
consistent with a Gyr duty cycle. 

Star formation in the context of bi-polar AGN feed¬ 
back triggered by cold gas was first tracked in |Li et al. 
(2015a), and that simulation exhibits episodic outbursts 
of star formation on Gyr time scales, as well as filamen¬ 
tary structures of cold star-forming gas tens of kpc long, 
while the minimum threshold of t coo i/tff varies from 5-20 
over the course of any single outburst. These simulations 
reproduce UV morphologies similar to those of coo l-core 
BCGs in CLASH, shown in Don ahue et al. fl2015|), and 
of lower-redshift BCGs studied in Tremblay et al. "([2015 ). 
The earlier simulation work described here inspire d the 
analytic preci pitati on model framework des cribed in Voit 
et al.| (|2015a|) and |Voit & Donahue| (|2015|) , which com¬ 
pared the implications of a minimum cooling time to 
free fall time to the bimodal entropy profiles seen in X- 
ray observations of clus te rs from the ACCE PT database 
(Cavagnolo et al. 2008). Voit et al. (2015) extends this 
precipitation framework to lower-mass galaxies and de¬ 
rives interesting implications for the connections between 
galaxy mass, mass of the central black hole, star forma¬ 
tion efficiency, and chemical composition. 

Since jet-triggered precipitation ought to have a mor¬ 
phological relationship with the jet and a characteristic 
timescale set by the AGN duty cycle, we can look for 
evidence that may support or contradict this prediction 
by comparing the SED-derived stellar parameter maps in 
RXJ1532 to X-ray measurements of recent AGN activity. 
In the single-aperture SED of RXJ1532, the starburst 
lifetime is log A4 [M©] = 8.8 ± 0.5 logioyr, so the bulk of 
the starburst has a lifetime on the order of ^1 Gyr. How¬ 
ever, the southern filaments and the northwest-southeast 
filaments are 10 7 — 10 8 yr old. The Chandra X-ray image 
reveals two well-defined cavities to the east and west of 
the BCG, and possibl e evid ence o f ghost cavities to the 
north and south (Hlavacek-Larrondo et al. 2013). The 
cavity refill times, whic h are the largest estimates of th e 
cavity ages provided by Hlavacek-Larrondo et al. (2013), 
are 6.3T0.7 and 8.2T0.7 x 10 Y yr. The cavities appear 
to be young relative to the timescale of the starburst we 
recover in our analysis, but match the ages of the young 
filaments. Our results in RXJ1532 are consistent with 
an ongoing process of clumps and filaments precipitating 
out of the ICM when pushed out of equilibrium by jets. 
Application of the SED fitting techniques developed in 
this paper to other BCGs in the CLASH sample will de¬ 
termine if this narrative is consistent for all of the star 
forming BCGs. 

Since jet-triggered precipitation ought to have a mor¬ 
phological relationship with the jet and a characteristic 
timescale set by the AGN duty cycle, we can look for 
evidence that may support or contradict this prediction 
by comparing the SED-derived stellar parameter maps in 
RXJ1532 to X-ray measurements of recent AGN activity. 
In the single-aperture SED of RXJ1532, the starburst 
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lifetime is log AG [M 0 ] = 8.8 db 0.5 logioyr, so the bulk of 
the starburst has a lifetime on the order of ^1 Gyr. How¬ 
ever, the southern filaments and the northwest-southeast 
filaments are 10 7 —10 8 yr old. The Chandra X-ray image 
reveals two well-defined cavities to the east and west of 
the BCG, and possibl e evid en ce of ghost cavities to the 
north and south (Hlavacek-Larrondo et al. 2013). The 
cavity refill times, whic h are the largest estimates of t he 
cavity ages provided by Hlavacek-Larrondo et al. (2013), 
are 6.3T0.7 and 8.2T0.7 x 10 Y yr. The cavities appear 
to be young relative to the timescale of the starburst we 
recover in our analysis, but match the ages of the young 
filaments. Our results in RXJ1532 are consistent with 
an ongoing process of clumps and filaments precipitating 
out of the ICM when pushed out of equilibrium by jets. 
Application of the SED fitting techniques developed in 
this paper to other BCGs in the CLASH sample will de¬ 
termine if this narrative is consistent for all of the star 
forming BCGs. 

The cavities in RXJ1532 also appear to be aligned 
with the northwest-southeast oriented young filament, 
and anti-aligned with the shape of the BCG filament net¬ 
work more generally. This morphological relationship is 
depicted in Figure [T6| The young filament traces one of 
the bright Ha filaments as well. The most prominent X- 
ray cavity corresponds to the brighter (western) end of 
this filament, suggesting a system with the western edge 
inclined towards us. Given the available data, we do 
not rule out a coincidence; however, the corresponding 
ages and morphologies suggest jet-triggered formation of 
the young filaments. We hypothesize that the northwest- 
southeast filament may have been the result of positive 
feedback triggered by compression from a jet inflating 
the X-ray cavities. 
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Figure 16. Left: The burst log lifetime map for the RXJ1532 
BCG, which emphasizes the young filaments, is shown with 
contours of 0.5-7.0 keV counts from the Chandra observation 
of the cluster shown in black. The red line traces the position 
of the young filament extending along a northwest-southeast 
axis. Right: The 0.5-7.0 keV Chandra image of RXJ1532. 
The red line is identical to the line shown depicted in the left 
panel. The red line points directly into the prominent cav¬ 
ity in the northwest of the image, and appears to lie along 
the axis connecting the northwest cavity to the much fainter 
southeast cavity, whic h is ev ident in the u nsharp-masked im¬ 


age presented in 


Hlavacek-Larrondo et al. 


(2013). 


The narrative we propose for RXJ1532 may be typi¬ 
cal for other BCGs in cool core clusters. The SFH of 
RXJ1532 agrees with a study of lin e em i tting B CGs in 
the SDSS survey produced by Liu et al. (2012). While 
their interpretation of the burst history differs from ours 
(they assume a stellar population divided into three 


components- a recent starburst, young stars, and old 
stars), they find that the majority of the flux they ob¬ 
serve in that sample of BCGs is due to stars forming 
within ^2.5 Gyr. Furthermore, estimates of molecular 
gas masses imply that B CG sta rbnrsts typically have 
fuel to last Gyr (|0’Dea et al. 2008]). While the case 
of RXJ1532 may be extreme in terms of SFR, its burst 
history may be the norm for star forming BCGs. The 
sum of the results descri bed here, taken i n ad dition to 
the evidence provided by Donahue et al. (2015), shows 
agreement between the resolved star structures we ob¬ 
serve in the CLASH BCGs and recent theoretical work 
on feedback in cool core clusters. 

If gas is condensing out of the ICM and feeding star 
formation, then we expect a general trend between the 
SFR and our M g . The idea of an AGN jet-driven mech¬ 
anism for ICM condensation could also account for our 
finding that as the SFR increases in the BCG, the SFR 
accounts for a larger fraction of the total cooling im¬ 
plied by M^. M^ r35 and M ^ ?t 20 can be interpreted as the 
cooling rate of gas in these radii in the absence of reheat¬ 
ing. Based on this interpretation, we suspect that larger 
starbursts occur in BCGs where gas is cooling more ef¬ 
ficiently. The correlation between cool core BCG SFRs 
and ‘efficiency’ can be explained neatly if cooling is lo¬ 
calized around AGN jets and cavities, or if in systems 
with stronger cooling a larger fraction of the gas with 
tcooi/tff < 20 is at t coo i/tff < 10. If larger SFRs oc¬ 
cur in BCGs exhibiting AGN feedback in larger areas, a 
larger fraction of unstable ICM gas inside a given radius 
may be triggered into condensing by feedback. This sce¬ 
nario is consistent with our finding that SFR is correlated 
with the area of the star forming re gions i n B CGs. 

In the feedback cycle modeled in |Li et al.| (|2015b|), a 
large amount of gas is cooled quickly with the initial on¬ 
set of jet feedback (on the order of 10s of Myr), and is 
slowly consumed by star formation (on the order of a 
Gyr). Larger SFRs tend to occur earlier in the feedback 
cycle, when proportionally more of the ICM is precipitat¬ 
ing. If this concern dominates the SFR-M g relationship, 
then we expect to see evidence that larger BCG star- 
bursts are younger on average. In RXJ1532, we see evi¬ 
dence connecting the durations of star forming knots and 
filaments to the duration of AGN activity (as evidenced 
by estimates of the ages of X-ray cavities in this system), 
although a more detailed study of multiple systems will 
be necessary to establish whether or not there is a rela¬ 
tionship between burst age, SFR, and low-entropy ICM 
gas. 


6. CONCLUSION 

We have conducted a detailed analysis of the star form¬ 
ing structure in the BCGs in the CLASH X-ray selected 
sample of galaxy clusters. Using the rich set of CLASH 
photometry, we estimated the dust reddening in BCGs 
with significant UV emission and calculated reddening 
corrected mean UV luminoisities and Ha+[N II] lumi¬ 
nosities. We compared these measurements to observa¬ 
tions of [O II] and H (3 taken using the Goodman spectro¬ 
graph. Additionally, we compared the UV derived SFR 
to X-ray properties calculated using the ACCEPT cata¬ 
log, including the core entropy, Kq, and predicted cool¬ 
ing rates for low entropy gas inside r = 35 kpc (M^ r 35 ) 
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and inside radii where t coo i/tff — 20 (M^^o)- We con¬ 
cluded our analysis by creating a resolved map of the 
starburst in RXJ1532.9+3021, for which we al so have an 
SPSS spectrum and d etailed X-ray data from |Hlavacek- 


Larrondo et al. (2013). 

Using measurements of [O II], [O III], and H/? lines 
in conjunction with broadband H<a+[N II] estimates, we 
constructed diagnostic diagrams for the CLASH BCGs 
in order to constrain the line-emission power source in 
these galaxies. Line emission in CLASH BCGs are pow¬ 
ered by a combination of star formation and a LINER- 
like source (possibly the signature of hot, young stars 
or interaction between the ICM and nebular gas), while 
the biggest starburst (MACSI931) has a line emission 
spectrum dominated by ongoing star formation. 

CLASH SFRs span a range of magnitudes up to > 100 
Mq yr -1 , and significant, extended star formation oc¬ 
curs in 10 out of 20 BCGs in our sample. Based on com¬ 
parisons with Ko and M^, we establish a link between 
the star formation in the BCG and the state of the sur¬ 
rounding ICM. All of the star forming BCGs with an 
SFR > 10 Mq yr 1 are consistent with a ^ 30 KeV cm 2 
entropy threshold, and a trend exists between SFR and 
M^. These findings imply SFR is fueled by a reservoir of 
low entropy gas. 

SEP analysis of RXJ1532.9+3021 reveals a long-lived 


starburst, with a log lifetime of 8 . 8 + 0 .5 logioyr, and a 
total star formation rate of II 8 I 42 5 M© yr -1 , which is 
consistent with our estimates from UV and Ra luminosi¬ 
ties. The overall burst timescale is much longer than the 
AGN on-cycle as inferred by the ages of AGN cavities in 
the ICM of this cluster, although several of the individual 
filaments are consistent with the ^ 60 — 80 Myr cavity 
refill times. These results are consistent with recent jet- 
triggered filaments super-imposed on an older long-lived 
starburst, which may have been the result of jets from 
previous AGN on-cycles. The burst history in RXJ1532 
is also c onsisten t with another study of stellar popula¬ 
tions by Liu et al. ( 2012 ) conducted on SPSS BCGs, 
so we hypothesize that m upcoming work we will find 
similar evidence for sporadic starbursts corresponding to 
episodes of AGN activity. 

The SFH of knots and filaments in RXJ 1532 suggest 
a jet- induced precipitation scenario such as Li et al. 
(2015b) is responsible for converting the ICM into cold, 
starforming gas. If true for all CLASH BCGs, this mech¬ 
anism would explain the relationship between the ther¬ 
modynamic state of the ICM surrounding CLASH BCGs 
and the SFRs in the BCGs. The increasing ‘efficiency’ of 
BCG SFRs relative to the cooling rates implied by M^ as 
a function of SFR is plausibly explained by this scenario, 
as well. 


APPENDIX 


This appendix includes the details of the parameter space chosen to fit the photometry of RXJ1532.9+3021 to a 
distribution of model SEPs. We describe the parameterization of the SFH, along with the parameter space we defined. 
Finally, we describe the consistency tests performed on fitting the CLASH SEP 


For the 16 band SEPs we const ruct using all the available bands of CLASH photometry, we used the Salpeter (1955) 
IMF and Bruzual & Chariot (2003) SSP. The SFH we fit consists of a uniform starburst imposed on a background 
population with an exponentially decaying SFR, thus the SFH is modeled by 


Ipe (t) = 


_ -Nearly -t/r 


e — (*age - a *{,)/t 


_ f FbM ea rly At 6 

I 0 


1p b (t) = 

'•Pnet ( t ) = 1p e (t) + Ipb (t) 


Gge t ^ Atfr 
Gge t + At^ 


(1) 

( 2 ) 

(3) 


which is a variant of the SFH described in 
population of stars, which is parameterized 


Moustakas et al. (2013). is the SFH for the background early-type 
3 y the time constant r. ^ is the burst SFH, and it is parameterized by 


the burst lifetime AG and the fractional burst amplitude F^. Our SFH consists of one burst for simplicity. 

We allowed the age of the galaxy to vary between 6 and 9.5 Gyr, and we allowed r to vary between | and ^ the 
age of the BCG, thereby ensuring that the background population corresponds to a quiescent, early-type galaxy. We 
allowed the metallicity of the stellar population to vary between 0.04 Z© and 1.6 Z 0 , and the dust attenuation Ay to 
vary between 0 and 5 mag. In order to sample a wide range of possible burst histories, we sampled burst parameters 
logarithmically, selecting A E in the range -3 < logA^ [Gyr] < 0.8, and F & in the range -2 < logF^ < 1.0. This 
parameter space is summarized in Table Al. We drew 10 4 models from this parameter space for each SED we fit to. 

SED Fitting Consistency Checks 

The equivalent width of [N Il]+Hcq EW([N Il]+Ha), is exqu isitely sensitive to the SF H of a galaxy, since it is a 
measure of the ratio of an SFR indicator to the red continuum (|Kennicutt||1998| |Leitherer 2005|). Since this value was 
measured directly by SDSS, and our model SEDs predict line strengths, we can compare the best fit EW([N Il]+H<a) 
for the CLASH SED with the SDSS measured value. This in turn indicates how reliable our estimate of the burst 
duration is. 

We perform this comparison by fitting the SED of fluxes extracted in a 3" diameter aperture centered on the 
coordinates of the SDSS fiber used to take spectra of RXJ1532. Our predicted EW([N Il]+H<a) = 269+120 A rest , 
which matches the SDSS measurement of 269.3+2.4 A rest . The probability distribution we recover also shows that 
the burst duration (AG) is logAG = 8.8 ± 0.5 [yr]. and burst mass (M&) is log M& = 10.83 ± 0.35 [M©]. The overlap 
between the SDSS fiber and the observable UV structure in RXJ1532 is substantial, so it is not surprising that we 
recover the same burst history and burst mass for this SED as the SED fit in the Results. 




























Table A1 

SED Fitting Parameter Choices 
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Parameter 

Range 

Units 

SSP 

BC03“ 


IMF 

Salpeter 6 


Model Draws 

10 4 


tage 

[6, 9.5] 

Gyr 

T 

[0.05, 0.2] 

tage 

Metallicity 

[0.04, 1.6] 

z© 

Ay 

[0.0, 5.0] 

mag 

log [O III]/H/3 

[-0.5, 0.5] 

dex 

log A t h 

[-3.0, 0.8] 

log Gyr 

log F b 

[-2.0, 1.0] 

dex 


a Bruzual fe Chariot] ( |2003| ) 

b Salpeter] ( |1955| > 



log EW([NII]+H<*) [A] 


Figure Al. We plot the distribution in A t b and EW([N Il]+H<a for all the models gridded in the parameter space summarized 
in Table Al. Each point corresponds to the values of these two parameters for each individual model, and the vertical line 
denotes 269.3 A res t, the SDSS measured EW. This plot reveals a densely populated trend tracing out a curve of decreasing 
burst lifetimes as a function of EW, with a wider, more sparsely populated envelope. 


As shown in Figure |AT| A t b is highly dependent upon the equivalent width. The sensitivity of the equivalent width 
to the burst history makes it straightforward to constrain a relatively narrow range of burst durations for the model 
SED. The agreement between the spectral equivalent width and our SED fit value is important because it shows that 
the H<a-f[N II] feature is detected strongly enough in the CLASH photometry that a meaningful constraint on the 
burst history can be made with it. 

We used iSEDf it to fit only the three bands of WFC3/UVISIS photometry that were used to estimate L^yy, and 
show that the results agree with our broadband estimates. We fit photomet ry extracte d from the region used to 
calculate L uv and assumed fit parameter constraints that are consisent with Kennicutt (1998). For this step, we 
assumed the metallicity in the stellar population is 1 Z 0 , that the observed age of the BCG, £ age , is at least 6 Gyr, 
and that the de cay ti mescale for the SFH , r, is 100 Gyr in order to force a continuous SFR model. We also assume a 


Salpeter (1955) IMF, Bruzual & Chariot (2003) SSP, and a Calzetti et al. (2000) reddening law. 

Using this method, we calculated an SFR of 99T24 M 0 yr -i , and an average reddening of E(B-V) = 0.26 db 0.04. 
As before, the SFR PDF is log-normal and we report the mode of the distribution. This result is consistent with the 
UV photometry results. The agreement implies that while taking a single-aperture SED fit washes out the correlation 
seen between the spatial distribution of the dust and the SFR, this effect is not a strong source of systematic variation 
in our fits. 
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